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ABSTRACT 



o 

^' 

o 

, The difference in the adiabatic indices of photons and non-relativistic baryonic matter in the early 

Universe causes the electron temperature to be slightly lower than the radiation temperature. 
Comptonization of photons with colder electrons results in the transfer of energy from photons to 
Ijr^ ■ electrons and ions owing to the recoil effect (spontaneous and induced). Thermalization of pho- 

tons with a colder plasma results in the accumulation of photons in the Rayleigh- Jeans tail, aided 

■ by stimulated recoil, while the higher frequency spectrum tries to approach Planck spectrum at 

■ the electron temperature T^"'*' = < T'™'"^'; i.e., Bose-Einstein condensation of photons occurs. 
We find new solutions of the Kompaneets equation describing this effect. No actual condensate 
is, in reality, possible since the process is very slow and photons drifting to low frequencies 
are efficiently absorbed by bremsstrahlung and double Compton processes. The spectral distor- 
tions created by Bose-Einstein condensation of photons are within an order of magnitude (for the 

■ present range of allowed cosmological parameters), with exactly the same spectrum but opposite 



in sign, of those created by diffusion damping of the acoustic waves on small scales correspond- 
ing to comoving wavenumbers 45 < < 10* Mpc"'. The initial perturbations on these scales are 
completely unobservable today due to their being erased completely by Silk damping. There is 
partial cancellation of these two distortions, leading to suppression of /j distortions expected in 
the standard model of cosmology. The net distortion depends on the scalar power index and its 
running dn^/dlnk. and may vanish for special values of parameters, for example, for a running 
spectrum with, n^, = l,dns/dlnk = -0.038. We arrive at an intriguing conclusion: even a null 
result, non-detection of //-type distortion at a sensitivity of 10 gives a quantitative measure of 
the primordial small-scale power spectrum. 

Key words, cosmic background radiation — cosmology:theory — early universe — 
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1. Introduction 



In the early Universe we have two strongly interacting fluids: cosmic microwave background 
(CMB) radiation and plasma (=ions+electrons). Electrons and CMB are strongly coupled due to 
the Compton process with the rate of energy transfer from CMB to electrons given by ^compton ~ 
1.4 X 10"^°(1 + z)'*s"'. This is much faster than the expansion rate , H{z) ~ 2.1 x 10"-"(1 + z)^s"' 
0, and keeps Ty - ^ Ty at redshifts z > 200, where Ty is the eff'ective temperature of radiation 
(to be defined below) and is the electron temperature. Ion s and electrons exchange en ergy via 
Coulomb colUsions at a rate /^coulomb ~ 3 x 10^^(1 + z)"'^^s ' ( Lifshitz & Pitaevskill981 ) keeping 



Te - Ti Tf., where is the temperature of ions. Thus we have Ty k x Ti to high precision at 
z > 200. Non-relativistic matter has an adiabatic index of 5/3, while the adiabatic index of radiation 
is 4/3. Adiabatic cooling due to the expansion of the Universe thus makes the matter cool faster 
than the radiation, while Comptonization tries to maintain matter and radiation at the same temper- 
ature by transferring energy from radiation to matter. As a result, the energy of radiation decreases, 
while the number of photons is conserved (negle cting absorption of pho tons by bremsstrah lung 



Zeldovich et al 



d 19681) and 



Peeblesl(ll968h and 



and double Compton). This eff'ect was discussed byi 
leads to a decrease in the temperature of plasma compared to the radiation temperature at z < 200, 
when Comptonization becomes inefficient because of depletion of electrons from recombination. It 
also drives the 21 -cm line spin temperature to below the radiation temperature, raising the possibil- 
ity of 21 -cm absorption of CMB photons by neutral hydrogen. This transfer of energy from CMB 
to matter happens at all redshifts, as long as there is significant free el ectron number density, and 



Chluba & SunvaevI (l2012l) 



results in deviations in CMB spectrum from the blackbody. Recently, 
have numerically studied these deviations from the blackbody in the CMB spectrum. Below we 
demonstrate that the problem of spectral deviations of CMB due to loss of energy to plasma, in the 
case of small spectral deviations and small Compton parameter y, has a simple analytic solution, 
when photon production and destruction is neglected. The resulting flow of photons towards low 
frequencies as the spectrum tries to approach the Planck spectrum is in fact the Bose-Einstein con- 
densation of photons. The Bose-Einstein condensation of photons unfortunately does not progress 
very far in reality, because low-frequency photons are efficiently absorbed by bremsstrahlung and 
double Compton processes, and the Compton process freezes out as the electron density falls owing 
to the expansion of the Universe and recombination. 

Early papers about energy release in the early Universe were concerned with exotic sources 
such as annihilation of matter and antimatter, primordial turbulence, decay of new unstable parti- 
cles, unwinding of topolo gical defects like d omain walls, and cosmic strings. Experiments begin- 
ning wiffi COBE FIR AS (iFixsen et al.lll996h have bee n unable to find a ny significant distortions 
in ffie CMB from blackbody. New proposals like Pixie ( Kogut et alJbol 1 ) are demonstrating that a 
tremendous increase in the sensitivity is possible in the future experiments. Pixie is proposed to be 
1000 times more sensitive than COBE FIRAS. There is hope that this is still not the last word, and 
even higher sensitivity might become possible in the future. Under these circumstances we decided 
to check what minim um levels of deviation s from blackbody spectrum are expected in the standard 



cosmological model (IKomatsu et alj|201 ih . In this paper we only consider the spectral distortions 



' We give the formula that is valid during radiation domination just for illustration. The Hubble rate remains 
much lower than the Comptonization rate at z > 200 and than the Coulomb rate at even lower redshifts. 
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arising before the end of the dark ages and beginning of reionization. We demonstrate that, in the 
absence of decay and annihilation of new unknown particles or any other new physics beyond the 
standard model, there are only three key reasons for significant global spectral distortions in the 
CMB. 

1. Bose-Einstein condensation of CMB photons due to the difference in the adiabatic indices of 
non-relativistic plasma and photons. 

2. The en ergy releas e due to dissipation of sound waves and initial perturbations due to Silk 
damping ( ISilklll968[) and the second-order Doppler effect due to non-zero pecuhar velocity of 
electrons and baryons in the CMB rest frame, both leading to the superposition of blackbodies 
( Zeldovich et al.|[l972 ) in the electron rest frame. 

3. The cosmological recomb ination radiation from hydrogen and hehum ( Dubrovichll975 



Peebles 



1968 



Zeldovich et al.ll 19681) 



We omit the distortions caused by the co smological recombination process (e.g., see 



Chluba & Sunvaevl[2006 : 



Sunvaev & Chlubal[2009l for more details and references therein) in the 



discussion below. They have narrow features and different continuum spectrum and can be distin- 
guished from y and ju distortions from the first two mechanisms. Cosmological recombination radi- 
ation can, however, become comparable to or larger than those discussed here dChluba & Sunvaev 



2012h . if the distortions from the first and second processes partially cancel. We also note that for 



any experiment with finite beam size there will be a mixing of blackbodies in the beam owing to the 
temperature fluctuations on the last scattering surface on scales sma fler than the beam size, lea ding 
to inevitable}; distortions of magnitude ~ (AT/Tf ~ 10"^ - 10 '° ( Chluba & Sunvaevll2004 ). 



In this paper we find new solutions for the Kompaneets equation describing the CMB spec- 
tral distortions arising from process 1 described above. A second mechanism results in a y- 
distortion and heating of the electrons at low redshifts. At high redshifts Co mptonization of CMB 



on the hotter electrons converts the y-t ype distortion to a /x-type distortion (iSunvaev & Zeldovich 



(ll970bl)Ullarionov & SunvaevI ( Il975bl) ). We demonstrate below that these two sources of distor- 
tions work against each other, individually resulting in spectral distortions with opposite signs but 
the same spectral shape in any part of the CMB spectrum. 

The cooling of electrons and the corresponding spectral distortions are easy to calculate and 
only depend on the standard cosmological parameters, such as baryon to photon number density 
nfi/wy and helium fraction, which decide the amount of energy losses by the CMB, the Hubble 
constant (Hq), and densities of constituents of the Universe, which in turn determine the expan- 
sion rate and thus the efficiency of Comptonization. On the other hand, the energy released by 
dissipation of sound waves crucially depends on the power in small-scale fluctuations, and thus 
the spectral index («s) in the standard cosmologi cal model, in addition to the other well-measured 
parameters of the standard cosmological model, ( Komatsu et al.ll20n ). It is interesting to note that 
for a spectrum with constant the energy release due to dissipation exceeds the energy losses to 
adiabatic cooling of baryons, and there is net heating of electrons. For a primordial spectrum with 
a running spectral index, the role of Comptonization with colder electrons and Bose-Einstein con- 
densation can become dominant, with the spectral distortions changing sign, and there can be net 
cooling of electrons and a corresponding decrease of CMB entropy per baryon (specific entropy). 
Additional increase in entropy during the recombination epoch due to superposition of blackbodies. 
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free streaming as well as Silk damping and the second-order Doppler effect only produce y-type 
distortions and can be distinguished from the ji distortions created in the earlier epoch. 

2. Thermodynamic equilibrium in the early Universe 

We can obtain the energy losses due to adiabatic cooling of baryons under the assumption that 
Compton scattering, bremsstrahlung, and double Compton scattering can maintain full thermody- 
namic equilibrium between the electrons and photons. This is true to high accuracy at z > 10^ while 
at lower redshifts, although bremsstrahlung and double Compton scattering cannot help in restor- 
ing a Planck spectrum at all frequencies, Comptonization keeps the electron temperature equal to 
the effective radiation temperature (cf. Eq. ( fTSI l) to high precision until z ~ 200. In the thermody- 
namic equilibrium, entropy is conserved, and we can use this to calculate the energy losses due to 
the expansion of the Universe without referring to Comptonization, since it is not really important 
which physical process is responsible for the equilibrium. 
We start with the thermodynamic relation, 

TydS = df/ + PdV - MdN, (1) 

where S is the entropy, U total thermal energy, P pressure, V volume, M chemical potential, and 
the number of particles. Ty is the common temperature of photons, ions, and electrons. We can 
ignore the last term for ions and electrons since their number is fixed after big bang nucleosynthesis 
end j^. For photons as well we can ignore the last term if we assume that their chemical potential is 
0, which is true in full thermodynamic equihbrium. With this assumption we can write the equation 
for total entropy per baryon, cr = i/«B, where s is the entropy density and ng is the baryon number 
density: 

, d(£/«B) + Pd(l/nB) 

do- = (2) 

-'r 



3 
2 

1 



E = a^T^ + -NnsksTy (3) 



P ^ -aj^T^y + NnsksTy. (4) 

Here or is the radiation constant and A' the number of non-relativistic particles per baryon, so that 
N - rif-b. We now integrate Eq. (|2|l to obtain 



-r^ +NkB In 

3«B 



hbC 



(5) 



where C is an arbitrary constant of integration, which is not important for our calculation. The 
first term above is the contribution from photons, and it is clear that in the absence of second 
term, Ty oc (l + z). The second term is the contribution of non-relativistic particles and causes 
the temperature to drop slightly faster in order to conserve the entropy per baryon. We can write 
Ty = T{1 + f), where T - 7'i( 1 + z)/( 1 -H z,) is the background temperature proportional to 1 -1- z, and 
f <*: 1 is the fractional deviation from this law. We can take the initial deviation f(z,) = at initial 



^ During recombination, the number of particles changes and the departure from thermodynamic equilib- 
rium also becomes significant. The following calculations are therefore strictly valid only before recombina- 
tion. 
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redshift Zi and then calculate the subsequent energy losses to adiabatic cooling of baryons at later 
redshifts. Thus we can write cr, correct to leading order in t, as 

4flRr3 4aKT\3t) , /r3/2\ 3NkBt 

0-=— + — +A^^Bln — +— r— . (6) 

3«B 3«B \«bC/ 2 

Equating the above to the initial entropy cr, = cr(z,) with t = 0,T = Ti, we get upon solving for t 
(and ignoring the last term in Eq. (|6]l compared to the second term) 

3A^«b^b , / 1 + z 
t T- In 



8flRr3 \l+zi 

^inflli 

1 1 + Zi 

5.9 X 10-1° . (l+z 



with AE/E - At, since Ey oc T'^, and Eb - jkBiif-bTg. we note that t is negative since z,- > z and 
there is net energy loss. \n{zmRx/Zmin) ~ 10 for the redshift range of interest, 200 < z < 2 x 10'', 
where a mixture of jj and y-type distortion s are produced. This result for the total energy extraction 
is consistent with the estimate obtained bv lChlubal (120051) and lChluba & SunvaevI (l2012h . 

To estimate the magnitude of total distortions, fi and y-type, we take Zmax = Z; = 2 x 10^, 
since at greater redshifts spectral distortions are rapidly destroyed by bremsstrahlung and double 
Compton processes and Comptonization. The latter redistributes photons over the whole spectrum. 
Forzmin - 200 wc get Ais/Zi - 5.4x10"^. This transfer of energy from radiation to baryons due to a 
difference in their adiabatic indices results in an inevitable distortion and Bose-Einstein condensa- 
tion of photons. We show below that the distortion has a magnitude of Fbec ~ (1/4) AE/E ~ 10~'. 
We define parameter Fbec and prove this result below. The corresponding deviation of the electron 
temperature from the equilibrium temperature of electrons in the radiation field is described by 
Eq.([T6]l and is {Ty - T^)/Ty ^ lO"'- at z = lO'' growing to ~ 0.01 at z = 500. 

In the Universe there are also inevitable processes like the dissipation of acoustic waves and SZ 
effect from reionization, which would lead to a normal SZ effect of similar or greater magnitude. 
Due to the additivity of small spectral distortions, these effects will cancel the Fbec distortions, and 
this diminishes the hope that this Bose-Einstein condensation will ever be observed. Nevertheless 
it is important to stress that in the standard cosmology there was an epoch lO*" > z > 1000 when 
Bose-Einstein condensation of CMB photons due to the difference in the adiabatic indices of matter 
and radiation was able to create a peculiar deviation of CMB spectrum from the blackbody. 



3. Alternative direct calculation of energy losses in CMB 

We would Uke to remind the reader that adiabatic cooling of baryons leads to the cooling of radia- 
tion, but Comptonization itself conserves the number of photons. We can write the equation for the 
evolution of average thermal energy density ^b in baryons before recombination as 

dEe _3 d«eb ^3 d^e 
dz 2 dz 2 dz 

The derivative in the first term on the right-hand side, d«eb/dz - 3«eb/(l + z), is just the decrease 
in number density of particles (ions and electrons) with redshift z due to the Hubble expansion, 
with «eb - Pb/jUmoi, whcrc pb is the baryon mass density and jUmoi the mean molecular weight. The 
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change in temperature has contributions from adiabatic coohng and also from energy gained from 
CMB by Comptonization: 

dre_^_^_ 

, ~ , r-ii ^ Compton \^) 

Az I + Z JATBHeb 

where 5compton - A-ksEynf^crTiTy - Ti^)lmscH{z){l + z) is the energy transfer rate per unit volume 
from radiation to baryons by Compton scattering, «e is the number density of free electrons, ctj the 
Thomson scattering cross section, k^, the Boltzmann constant, Ey the energy density of radiation, c 
the speed of light, me the mass of electron, and H{z) the Hubble parameter. That redshift decreases 
with increasing time, so the terms with "+" sign are cooling terms and terms with "-" sign are 
heating terms. The change in photon energy density can therefore be written as 

dEy _ I dEy \ / dEy \^°"""°" 

"dr~\"drj '^['dz) 

_ 4Ey 

— Z 1" S Comptom (10) 

1 + z 

where the first term is just the adiabatic cooling due to the expansion of the Universe. We can 
estimate the energy transfer, S'compton, by noting that the baryon temperature T^. = Ty cc I + z to 
a high accuracy until z ~ 200, where Ty = 2.725(1 + z) K is the CMB temperature. Using this to 
evaluate the total derivative of electron temperature dT^/dz on the left-hand side of Eq. ^ and the 
first term on the right-hand side, we get 

3 2Ty 3 dr. 



^Compton - -^BWeb-; — - -^B«eb-7^ 

2 1 + Z 2 dz 



(11) 



(12) 



Thus we have fractional rate at which energy is lost by radiation to baryons 

(d£y/dz)^°'"P'°" ^ 5 Compton ^ E^lEy ^ 5.9 X 10''" 

Ey Ey (1+z) 1+z ■ 

We can now calculate the total fractional energy losses of radiation that contribute to the spec 
tral distortions between redshifts Zmin and Zmax, 

dE. 



AE _ p'""" dEy 

~ Ey[dzj 



= -5.9xlO-'°ln(i^^V (13) 

\ 1 + Zmin / 

We have integrated dE/E, since immediate distortions are proportional to AE/E, and the distortions 
can be added linearly if they are small. This is the same result as Eq. Q. 



4. Kompaneets equation 



The interaction of radiation with electr ons through Compt on scattering or Comptonization is de- 
scribed by the Kompaneets equation (iKompaneetsI 1 1 9561) in the Fokker-Planck approximation, 
when the energy transfer in each scattering is small compared to temperature, and the incoming 
photon distribution is wide compared to the width of the scattering kernel: 



dn 1 d 
— = — — + «' + 
oy ox 



Tf. dn 
^BB dx 



(14) 
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where we have defined the Compton parameter 



y(Z, Zmax) 



dz- 



(15) 



which is convenient to use instead of time or redshift. We start our calculation at the reference 
redshift Zmax- One can furthermore introduce the Compton equilibri um electron temperature in a 



radiation field ( Levich & Sunvaevlll971 



Zel'Dovich & Levichlll97C ) 



qu 

% 



Tl'^ /«(l+n)x4dx 



(16) 



4 J nx^dx 

Here Tbb is a reference temperature which is equal to the radiation blackbody temperature if the 
initial radiation field is a Planck spectrum, and x = /!v/A:b7'bb is the dimensionless frequency with v 
the frequency of photons. The equations written in this way factor out the expansion of the Universe 
and are applicable to cosmology, as well as non-expanding astrophysical systems. This equation 
and the Kompaneets equation form a coupled system to be solved simultaneously. 

The three terms in the inner brackets of Eq. (fT4l i describe recoil (n), induced recoil («-), and 
the Doppler effect of the thermal motion of electrons (Tg/TsBdn/dx). For y <K 1 and an initial 
blackbody spectrum with temperature Tbb, »bb = l/(g' ~ 1). a s imple first-order correction to 
the blackbody spectrum can be found ( Zeldovich & Sunvaev 1969 ). By approximating^ « + n~ x 
-dn/dx in Eq. (fT4l i we get 

4<9n(jc,y) 



dn(x,y) 1 d 
Ady x^ dx 



dx 



(17) 



where A = T^/Tsb - 1- This particular form of equation emphasizes that the Doppler term 
dominates over recoil and gives rise to the spectral distortion by boosting low-energy pho- 
tons to higher frequericy. Th e above equation can be reduced to the diffusion or heat equation 
( Zeldovich & Sunvaevll969 ) by changing variables, for which the Green's function is well known. 

For small distortions, the solution takes a particularly simple form and can be found by substi- 
tuting ^BB on the right-hand side of Eq. (fT4l i or Eq. ( fTTl i: 

+ 1 



An = n{x, Ysz) - nBsix) = i'sz" 



1 



(18) 



where we have defined Ysz - Jq Ady. The SZ parameter is related to fractional energy release by 
the formula Ys,z - {^/4-)AEy/Ey. This can be obtained by integrating the above spectrum over all 
frequencies. 



5. New solution of Kompaneets equation for Comptonization of photons with 
colder plasma 

The Kompaneets equation describes, to lowest order, the competition between the Doppler boosting 
of low-energy photons to high energies and the down scatter of high-energy photons to low energies 
by recoil and stimulated recoil. The initial radiation spectrum is blackbody with temperature Tbb- 
If the electrons are also at the same temperature, Te = Tbb, then these two effects cancel each 
other out exactly. If > Tbb, the Doppler boosting becomes stronger than recoil and we have the 
normal SZ effect in the limit of small y parameter. If Tg < Tbb then Doppler boosting is weaker 

^ For a Bose-Einstein distribution with constant chemical potential this expression is an identity. 
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Fig. 1. Difference in intensity from the blackbody radiation for normal SZ effect witfi Fsz = 10 ^ 
and negative SZ effect with Fbec = 10"''. 



and we have a net movement of photons from high to low energies, which is the complete opposite 
of the SZ effect. But since at - T^b the two effects exactly balance each other out, at linear 
order the spectral distortions for < Tbb would be exactly the same as for Tg > Tbb, just with the 
opposite sign. 

For Tf. < Tbb, we can write the approximate Kompaneets equation by approximating dn/dx x 
-(n + r?) in Eq. (fT4l l. 

dn 

Analytic solutions can be obtained by the standard method of characteristics for the recoil 
(n) and induced recoil («^) terms individually, and they cons ist of photons moving down the fre 



(1 



Tbb x~ dx 



(19) 



1971 



lUarionov & Svunvaev 



quenc y axes at a speed proportional to x^ for the recoil term (lArons 

1972 ) and at a speed pr oportional to the x^n for the induced recoil term ( IZel'Dovich & Levich 



(ll969l)lSvunvaevl (1197 II) V To get the linear solution, we can just substitute n{x) - «bb W on the 
righthand side of Eq. ( fT9] l and define the parameter for the amplitude of distortion in the case 
Te < Tbb as Fbec - Jq'CI - 7^e/rBB)dy, the result is Eq. (fTST l with Ksz = -^^bec- We should empha- 
size that the linear solution, Eq. ( fTST i, can be arrived at in a straightforward way by just substituting 
n{x) - riBBix) on the righthand side of the full Kompaneets equation, Eq. (fT4l i. without any restric- 
tion on Tg. However the dominant physical effects when Tg > Tbb (Doppler) and when < Tbb 
(recoil and induced recoil) are completely different and that in equilibrium they balance each other 
out exactly gives us the same mathematical form of the linear solution in both cases. 

In the Rayleigh-Jeans part of the spectrum we now have an increase in the brightness tempera- 
ture of NT IT - +2yBEC, which is independent of frequency and thus maintains the Rayleigh-Jeans 
shape of the spectrum. In the Wien part we have a decrease in intensity A/// - An/n^ - x^Ybec- 
Figure [U shows the difference in intensity from the blackbody radiation for normal SZ effect with 
Ysz = 10"^ and the negative SZ effect with Fbec - 10"^. Figure |2] shows the spectrum that would 
be achieved at a high value of the Compton parameter y. It is a Bose-Einstein spectrum with the 
dimensionless chemical potential^yu defined by the equation with «(v) = i/(gM':B7',,+/j _ i)^ where 



This definiti on has a sign difference with respect to the usual definition of the chemical potential in 
thermodynamics i Landau & Lifshitzll980 ) used in Eq. l[T). 
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Fig. 2. Difference in intensity from the blackbody radiation for distortion defined by the equation 
with n(v) = \ i^ei^^lk-eJ^+ii _ the temperature to which the spectrum relaxes for an initial 

energy addition/loss ofAE/E with AT^JT = 0.64A£'/£' = +2.5 x lO"''. AE/E = ±4 x 10"'' is the 
energy addition/loss that gives rise to SZ distortion of 10"^. The dimensionless chemical potential 
fi is given by ju = 2.2AT/T = ±5.6 x lO"'*. 



Tfj is the temperature to which the spe ctrum relaxes for an initial en ergy addition/loss of AE/E with 



AT^/T = QMAE/E = ±2.5 x 10 '' (llllarionov & Sunvaevl[i975ah . and AE/E = ±4 x 10 is the 



energy addition/loss that gives rise to SZ distortion of 10"^. The chemical potential is given bjO 
fi = 2.2AT/T - ±5.6 x lO"**. This is the spectrum that an initial spectrum with Ibec, i'sz = 10"^ 
will approach at high y at ~ 10 The frequency at which the distortion crosses zero is at 
jc = 2.19 compared to x = 3.83 for the SZ distortion in Fig.[Tl 

The y-type and yU-type distortions expected in the early Universe ( calculated in the later sec- 



2006) 



tions) are compared with the cosmological recombination spectrum (IRubino-Martm et al. 
in Fig. [3] Clearly ji type d istortions have a different spectral sha pe than the recombination r adia 
tion (both from hydrogen Jchluba & Sunvaevll2006l) and helium (iRubino-Martm et al 



20081) ') and 



y-type distortions and can be distinguished from the last two. This is very important because the 
information in the fi type distortions about the early Universe physics can be extracted. On the other 
hand, the y-type distortions from the early Universe get swamped by the much larger y-type dis- 
tortions from the low redshifts and the two contributions are difficult to separate. The fi distortions 
expected from the early Universe also have higher magnitude than the recombination spectrum in 
the Rayleigh- Jeans part of the spectrum but has no quasi-periodic structure like the cosmological 
recombination radiation. 

This flow of photons toward lower frequencies as the spectrum tries to approach the 
Planck spectrum due to recoil, and induced recoil is Bose-Einstein condensation of photons 



(llllarionov & Sunvaev 



1975a : 



Landau & Lifshitz 



19801) . We show the evolution of the spectrum 



(solution of the coupled system of Eqs. (fT4l l and ( fT6b ) starting with the initial distortion given by 
Eq. ( fT8b with -Fsz = isEC = 10"^ in Fig.|4] In the Rayleigh- Jeans region intensity ly oc T^s and 
therefore AT /T - AI/I. The initial evolution is similar to that of a spectrum with positive isz with 
the photon distribution approaching a Bose-Einstein spectrum defined by n(v) - l/(e'"'/^B7'c+/j _ j^-j 



' These relationships can be easily derived using photon number and energy conservation and requiring that 
the final spectrum have the equilibrium Bose-Einstein distribution. 
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Fig. 3. Comparison of positive and negative ^-type distortions expected in the early Universe (to 
be calculated in later sections) with the y-type distortions before reionization and the cosmol ogi- 
cal recombination spectrum (taken from 
illustration. 



Chluba & Sunvaev 



2006 : 



Rubifio-Martm et al 



2008) for 




1e-06 



1e-05 0.0001 0.001 0.01 

X 



100 



Fig. 4. Evolution of initial spectrum aty = given by Eq. ( fTSl l with -Ksz - ^bec = 10"^. \AT/T\ = 
\TeslT - 1 1 is plotted, where Teff is the temperature of a blackbody spectrum with the corresponding 
intensity at a particular frequency. Bose-Einstein spectrum defined by n{v) - l/(e'"'/'^B7'c+/; _ j^-j 
with negative /i is labeled '-ve BE', red dashed, double-dotted line in the figure. The distortions 
are positive in the Rayleigh-Jeans region (labeled '+ve') on the left side and negative in the Wien 
region (labeled '-ve') on the right side of the figure. 



with negative // (marked '-ve BE' in the figure). We should emphasize that there is no singularity in 
the actual solutions of the Kompaneets equation plotted above. The singularity is just in the math- 
ematical formula, which correctly describes the spectrum at high frequencies, x » The actual 
spectrum deviates from the Bose-Einstein spectrum near the singularity (positive everywhere in 
the Rayleigh-Jeans region) and can be described by a chemical potential decreasing in magnitude 
with decreasing frequency. The evolution at y > 1 is therefore very different from the positive 
Fsz case. For x » |ju| the spectrum is the Bose-Einstein spectrum at the electron temperature. For 
X «c |yu| there is an excess of photons compared to the Bose-Einstein spectrum which grows with 
time, a feature of Bose-Einstein condensation. The fractional change in temperature AT/T takes 
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1e+06 1e+07 
energy injection redshift (Zj) 

Fig. 5. Regions of the blackbody photosphere probed by different experiments, COBE FIRAS with 
sensitivity /i = 10"^, Pixie with sensitivity /i = 10"** and a hypothetical experiment with sensitivity 
of ju = 10"^. Shaded regions below the curves are allowed and inaccessible to these experiments. 
We also plot the ratio of final ju distortion today to the initial // distortion as a function of en- 
ergy injection redshift. We define the blackbody surface as the energy injection redshift such that 

^final /^initial ^ Ije^Zo, = 1.98 X 10^ 



a very simple form in the Rayleigh-Jeans region (x <*: 1) for a Bose-Einstein spectrum and can 
be understood as follows. It is given, in terms of frequency referred to the electron temperature, 

Xe = liv/kT,,, by 

— = ^- - 1 = (20) 

For |yu| «: .)c^ «; 1 we have AT/T - -fi/Xf,. Thus the fractional temperature deviation is positive 
for negative fi and negative for the positive fi. For x^, greater than and close to |//|,// < the spec- 
tral distortions can become very large and exceed unity. Since the high-frequency spectrum is in 
equilibrium at the electron temperature by y = 10, the subsequent evolution is very slow, and \fi\ 
decreases while the electron temperature increases slowly as the photon distribution approaches a 
Planck spectrum with the extra photons accumulating at low frequencies; i.e., Bose-Einstein con- 
densation happens. The low-frequency spectrum approaches the stationary solution of Kompaneets 
equati on with only the induced recoil term, n{x) oc l/x^, with a continuous flow of photons towards 
x^O Jsvunvaevll97ll) . Then, TJT- 1 just increases from -2.558 x 10"^ to -2.555 x 10"^ in going 
from y = 10 to = 100 for the chosen energy losses of AE/E = 4 x 10"^. The effect discussed 
below in the real Universe is of magnitude ~ 10"^. 



6. Bose-Einstein condensation of CIVIB photons in the early Universe. 



In the early Universe, before z ~ 2 x 10^ double Compton and Comptonization destroy any spec- 
tral distortions and maintain the Planck spectrum of CMB. For small distortions due to single, 
quasi-instantaneous episode of energy release, we can write the ratio of final-to-initial fi as an ex- 
ponential functi on of redshift defined by the square root of the product of Comptonization and 
absorption rates ( Sunvaev & Zeldovichll 1970b ). For a double Compton process, this formula gives 
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(iDanese & de Zottill 19821) 



, .final 

G(z) = L . (21) 

where Zdc ~ 1.98 x 10'' defines the "surface of the blackbody photosphere". We call this function 
G(z) the blackbody visibility function for spectral distortions. We call the region z > Zdc where the 
initial // distortion can be reduced by a factor of more than e the "blackbody photosphere". Thus 
inside the blackbody photosphere, a Planckian spectrum can be established efficiently. Figure |5] 
shows how the ju distortion at high redshifts decreases due to the double Compton and Compton 
scattering. The regions of the blackbody photosphere allowed by different experiments are shaded. 

At redshifts 10^ < z < 2 x 10^, we have the Compton parameter y > 1, and the Bose-Einstein 
spectrum with a negative chemical potential is established at x > 0.010. At lower frequencies, 
bremsstrahlung and double Compton create a Planck spectrum corresponding to the electron tem- 
perature. At z < 10"* the Compton parameter y < 0.01 and the distortions created can be de- 
scribed by the linear solution of Eq. ( fTSl ). Thus the net distortions created in the early Universe 
owing to adiabatic cooling of baryons are a linear superposition of distortions corresponding to 
different values of y, and a few of them are shown in Fig. |6l We show the evolution of actual 
spectrum rEq.[T4land[T6]l including the effects of bremsstrahlung and double Compton proce sses 



dBurigana et aljl991 



Chluba& Sunvaevl2012 : 



iijl 



Hu&Silkll993 



Illarionov & Sunvaevll975a|[bl) in 



the ACDM cosmology with WMAP parameters (iKomatsu et al 



201 11) . The reference spectrum is 



the blackbody CMB spectrum at z = 10^. At x < 0.01 the spectrum is dominated by bremsstrahlung 
and double Compton processes that destroy the condensing low-frequency photons. The effect of 
bremsstrahlung compared to Comptonization becomes stronger with decreasing redshift and can 
be felt at successively higher frequencies. 

At z > 5 X 10"* the spectral distortions in the x > 0.01 region can relax to a Bose-Einstein 
spectrum due to high value of parameter y. Contributions from lower redshifts are described by the 
analytic solution of Eq. ( fTSl l. The distortions created at z < 10"* are just the negative of the usual SZ 
effect and would be completely overwhelmed by the SZ effect from reionization at z ~ 10 which 
is expected to be Fsz - {kBTelmf.c^)Tn ~ 10"^ x 0.1 = 10"^, where Tn ~ 0.1 is the optical depth 
to the last scattering surface due to reionization and Tg ~ lO^K is the average electron temperature 
during reionization. The distortions created by energy loss at z > 10^ can, however, be destroyed 
only by energy injection at z > 10^' since at lower redshift Comptonization cannot create yu-type 
distortions. 



7. Energy release from dissipation of acoustic waves during radiation 
domination 

In the standard cosmology the main source of energy injection at high redsh ifts is the dissipation 
of acoustic waves in the baryon-photon pla sma (ISunvaev & Zeldovichll970bl) due to Silk damping 



(lSil3[l96i^ . 



Sunvaev & Zeldovich 



(Il970al) first proposed using the resulting spectral distortions to 



* Strictly speaking, a Bose-Einstein spectrum is established at y ~ few as can be seen from Fig|4]However 
at y ~ 0.25 - 1 the spectrum is already very close to Bose-Einstein at high frequencies. We use this to divide 



the energy release into //-type and y-type estimates. Exact results are presented elsewhere jChluba et al.l2012h 
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Fig. 6. Evolution of spectral distortions in CMB in ACDM cosmology. The values of total Compton 
parameter^ = y(100,z) for different curves ai-e (z - 1.5 x 10^, v = 109), (z - 4.6 x 10^, y - 10), 
(z= 1.5xlO\y= 1), (z = 5.9 xl0^y = 0.16), (z = 1.2 x 10^y = 0.005), (z = 1760,y = 3 x 10"^), 
(z = 840, y = 1.8 X 10"*^) and (z = 220, y = 1.2 x 10""'). Reference spectrum is the blackbody 
CMB spectrum at z = 10^. At low frequencies the spectrum is dominated by double Compton 
and bremsstrahlung which maintain Planck spectrum at the electron temperature. The zero point 
at low frequencies moves to the right at low redshifts because bremsstrahlung dominates over the 
Compton scattering at higher and higher frequencies and absorbs photons, although both are slower 
than the expansion rate. 



measure the spectral index and power in the small-scale fluctuations, which themselves do not sur- 
vive. Later during recombination the second-order Doppler effect due to non-zero electron/baryon 
velocity in the CMB rest frame also becomes important, but at that time only y-type distortions 
similar to the SZ effect arise, and jj distortions cannot be create d. We can estimate the amount 
of fi distortions created by dissipation of acoustic waves following lHu et al.l (1 19941) by calculating 



how the total power in the density fluctuations changes with time because of photon diffusion. The 
comoving energy density in acoustic waves in the photon-baryon plasma (neglecting the baryon 
energy density) is given hy Q - pyCj{(5y(x)^), where py is the comoving photon energy density, 
~ 1/3 is sound speed squared, 6y is the photon density perturbation at position x, and angular 
brackets denote the ensemble average: 



r d^k 

J (2^ 



(22) 



where Py(k) 



A^Xk)PL(k ), and P'(k) is the initial power spectrum. It was shown by 



Mukhanov & Chibisovl (119811) that primordial fluctuations from inflation can have a spectrum de- 
viating from the scale-invariant Harrison-Zeldo vich spectrum, with spectral index « , < 1 . In the 
cyclic ekpyrotic models > 1 is also possible (iLehners et al.ll2007h . WMAP gives the constraint 
on curvature perturbation in comoving gauge, ^. This is related to the gravitational perturbation i// 
in the radiation era (assuming neutrinos are free streaming) by the relation i// = ^/(2/5Rv + 1.5), 
where Ry = Pv/iPy + Pv) ~ 0.4, py is the neutrino energy density and 6'^ - -2i^ dMa & Bertschingei 



I995h . Thus P'y = 4/(2/5/?,, + l.5 f Pr ^ iA5Pr and P, = (A,2n^ /fYk/ k)) 



(is-l + ^Klnjl/A-o) 



0.002 Mpc-', Af = 2.4 x 10"'' ( dKomatsu et al 



2011 



R. Keisler et al. 



20111) ). r = dns/dlnk 



is the running of the index. The transfer function for modes well inside the horizon before recom- 
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bination is given by 



and the diffusion scale by (lKaiseilll983 



^ 3 cos(;tr,)e"*'/*°" 



Weinberg|l2008) 



d(lAD^) 



1 _ r 



dz 



c(l+z) 



16 

+ 



16 



6Hil + R)n^o-T \l + R 15 
c(l +z) 



(23) 

(24) 
(25) 



dz 6H{1 + R)n^o-T \l + R 15/' 

where R = 3ph/4-py, pi, is the baryon energy density. 

Replacing cos^{krs) with its average value over an oscillation of 1 /2, we get the energy release 
per unit redshift 



dQ/dz 
Py 



■y ^ dz 



-1 r dH 

3"J (2^ 
4.3Af d(l/^D') 



_2,2/,,2d(l/^D^) 



dz 



dk 



dz I 

4.3A^d(l/fcD')2_(3,„,)/2^^,,.ipJ 



n + 1 



dz 



(26) 



For a running index, the above integration must be done numerically. We are interested in radiation- 
dominated epoch where // distortions are generated, and we have H(z) = HoD.r"(l + zY and 
ndz) - («HO+2nHeo)(l+z)'' = «eo(l+z)"'- In addition we have /to = A^'^^(l+z)^/2 andd(l//tD^)/dz = 
-3Ad{1 + z)"^, and the above equation simplifies to 



dQ/dz 

Py 



UA. 



£2-(3+«J/2p 



+ 1 



4(l-«s)/2 



(1+z) 



(3ns-5)/2 



where 



An = 



8c 



135i/ofiy^neOO"T 



= 5.92 X 10"'Mpcl 



(27) 



Hu et al 



(1994). The diffusion 



The redshift and «j dependence obtained above matches that of 
scale /to = 46Mpc"' at z = 5 x 10^ and fco = lO'^Mpc"' at z = 2 x 10^ Thus we are probing the 
primordial fluctuations on very small scales that are not accessible in any other way. 
Thus for a Harrison-Zeldovich spectrum with n^-\ 

dg/dz -7.8 X 10-'' 



Py 

iZ=5xI0-' 
^=2x10" 



1 +Z 

2.9 X 10"' 



Py 



(28) 



Chluba & SunvaevI (120121) . For 



We note the surprising fact that this redshift dependence is exac tly the same as for the en ergy 
losses due to adiabatic cooling Eq. ( fTlT l, as already pointed out by 
the currently preferred value of - 0.96 we have 

de/dz -1 X 10-*^ 



Py 

^:=2xlQ<' 

Py 



(1+2)1.06 
1.8 X 10^1 



(29) 
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Figure |7] compares the rate of energy release multiplied by the blackbody visibility function, 
G(z)(l + z)dQ/dz/Ey for different initial power spectra with the energy losses due to the adiabatic 
coohng of baryons, (1 + z)S compton / Ey obtained by solving Eq. (|9]i for the standard recombination 
history dSeager et alj I200C ) calculated using the effective multilevel approach and taking recent 



corrections into account (lAli-Haimoud & Hiratall201 1 



Chluba & Thomaj201 11) . The spike due to 



increase in Silk damping during recombination is clearly noticeablqj. The smaller spikes due to 
Hell and Hel are also noticeable. Recombination of each species leads to a decrease in the number 
of particles, hence to adiabatic losses. But the photon diffusion length, and the associated energy 
release, increases owing to a decrease in the electron fraction. The transition from radiation to mat- 
ter domination results in unnoticeable change in the slope of the curves. Even after recombination 
photons continue to mix on horizon scales because of free streaming, while the energy losses from 
Comptonization drop sharply as a result of depletion of electrons. We switch to the free-streaming 
solution described in appendixlAlwhen diffusion length In/ku becomes equal to the comoving par- 
ticle horizon. In the blackbody photosphere region, the blackbody visibility function G{z) makes 
the curves drop sharply. Global energy release and y-type distortions during the free streaming 
epoch after recombination are due to the superposition of blackbody spectra on horizon scales. 
The energy injected into y-type distortion during recombination in the peak (800 < z < 1500) is 
AE/E ~ 10 for - 0.96 (see also footnote|7]). The y distortion due to mixing of blackbodies 
during free streaming after recombination and up to reionization , 20 < z < 800, calculated using 
the result in the appendix|Al is AE/E ~ 10"^. These are, however, difficult to separate from much 
the larger thermal and non-linear Doppler y distortion from reionization of Ysz ~ 10"^. 

Energy releases for different power spectra, including the ones with running spectral index, 
are summarized in Tables [1] and |2l The energy losses from adiabatic cooling of baryons during 
the same time period is ~ 2.2 x 10"^. This equals the energy release, for example, for a spectral 
index of «s = 1.0,dns/dlnk = -0.038 for a running spectrum. The net /j. distortion as a function 
of spectral index «s without running is shown in Fig.[8]and for «s - 1.0, 1.02, 1.05 as a function of 
the running of the index in Fig.|9l For low values of «s and high negative values of running, Bose- 
Einstein condensation dominates, and net distortion approaches a constant value of -3 x 10"^. A 
detection or non-detection at a sensitivity of 10"^ of ju-type distortion is thus a very sensitive probe 
of primordial power spectrum on small scales. Given the importance of the effective heating rate for 
the amplitude of the net distortion, it is impor tant to carry out a m ore careful and refined calculation 
of the problem, which we present elsewhere ( Chluba et al.ll2012 ). 

Over the years a number of ground- and balloon-based experiment s that aimed at mea- 
suring the absolute CMB brightness have been conducted before COBE ( Johnson & Wilkinson 



1987 



Levin et al 



1992) and after it ( Bersane 



Raghunathan & Subrahmnavan 



200C; 



Stages et al 



li et al 



1996c 



994 



Fixsen et al 



Zannoni et al 



2011 



2004allb 



i 



20081) . The best Hm- 

its on CMB u-type spectral distortions from these experi ments are at the level of 10 - 10"^ 



dFixsen et al 



1996 



Gervasi et al.ll2008 



Seiffert et al.ll201 Ih . Achieving a sensitivity of 10 would 



require subtracting foregrounds due to synchrotron emission, free free emission, dust e mission, and 



spinn ing dust emission at the same precision level. The proposed experiment Pixie jKogut et al 



201 11) aims to achieve this goal by using 400 effective channels between frequencies of 30GHz 



^ The actual energy injected in the recombination peak would be much less because of the transition to free 
streaming. 
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Fig. 7. Sketch of fractional rate of energy release due to Silk damping and free streaming for dif- 
ferent initial power spectra. Also shown for comparison is the rate of energy loss due to adiabatic 
cooling of baryonic matter. 
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Table 1. Energy injection in fi distortions during 5 x 10"* < z < 2 x 10^ for different initial power 
spectra without running compared with energy losses due to Bose-Einstein condensation. 




Fig. S. fi - 1 AAE/E as a function of spectral index rig without running. 



to 6THz. Their simulations indicate that an accuracy of InK in foreground subtraction is achiev- 
able. There is, however, some uncertainty in our understanding of the foregr ounds and possible 



systematics as indicated by an observed excess signal at 3GHz by ARCADE ( IFixsen et al 



201 ih 



experiment, which is not completely explained by the current galactic and extragalactic emission 
models; nevertheless achievement of a goal of 10"^ in sensitivity, although challenging, seems 
possible in the near future. 
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Table 2. Energy injection in distortions during 5xl0'*<z<2xl0^ for different initial power 
spectra with running compared with energy losses due to Bose-Einstein condensation. Energy in- 
jection values < BEC value in magnitude are shown in red. 
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Fig. 9. ji - \AlS.EIE as a function of running of the spectral index r for three values of spectral 
index n^. 




8. Conclusions 



Comptonization of CMB photons with slightly cooler electrons results in spectral distortions that 
then evolve (in the absence of any other physical process) towards the Bose-Einstein condensate so- 
lution, which is the equilibrium solution and consists of a blackbody spectrum with zero chemical 



potential, with the extra photons accumulating at zero frequencies (llUarionov & Sunvaev 



Landau & Lifshitz 



1975a : 



19801) . This is because the radiation has more photons than can be accommo- 



dated in a Bose-Einstein spectrum with non-negative chemical potential at the electron temperature. 
However, in reality, bremsstrahlung and double Compton scattering destroy extra photons at low 
frequencies. Thus in practice we have the high-frequency spectrum evolving towards a blackbody 
with the extra photons slowly moving down in frequency to be eaten up by bremsstrahlung and 
double Compton scattering. This is what is seen in the numerical solution in Fig.|6] 

The difficulty of conserving photon number is the reason Bose-Einstein condensation (or ac- 
cumulation of photons at low frequencies) is not observed in astrophysical systems. Bose-Einstein 
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Klaers et al 



conde nsation of photons has only been achieved in the laboratory very recently by 
(I2OIOI) . In the real Universe we also have inevitable sources of energy injection into the CMB, 
such as the SZ effect after reionization, which creates distortions of similar magnitude but with op- 
posite signs canceling and overwhelming this effect. The situation is much more interesting for // 
distortions, which can be created only at high redshifts and thus have no foregrounds at lower red- 
shifts. In standard cosmology, positive yU distortions are created by dissipation of acoustic waves on 
small scales due to photon diffusion. For the currently allowed values of cosmological parameters, 
there is a possibility that the positive ^ distortions from photon diffusion can almost exactly cancel 
the negative distortions from Bose-Einstein condensation of CMB, leading to a net distortion 
that is much smaller. Bose-Einstein condensation can even dominate, leading to a net distortion 

/Y 3 X 10 to which Silk damping contributes, but only as a small perturbation. Nevertheless, 

we must emphasize that the photon number conserving Comptonization, along with differences in 
adiabatic indices of radiation and matter, creates a unique system in the early Universe in which 
photons can begin to evolve towards a Bose-Einstein condensate. It is a remarkable coincidence 
that a completely unrelated physical process, diffusion damping of primordial perturbations, can 
produce distortions of almost exactly the same or greater magnitude, leading to the suppression 
of this unique effect in astrophysics. We thus arrive at an intriguing conclusion: even a null result, 
non-detection of ju-type distortion at a sensitivity of 10"^, rather than just placing a upper limit, ac- 
tually gives a quantitative measure of the primordial small-scale power spectrum. The importance 
of improving the experimental sensitivity to reach this critical value cannot be overemphasized. 



Appendix A: Erasure of perturbations due to free streaming 

We can derive the solution for erasure if perturbations on horizon scales due to free streaming as fol- 
lows. We start with the Boltzmann hierarchy, ignoring quadrupole and higher order moments. We 
define multipole(^) moments of the temperature anisotropyAr/r s = 2;('~'y(2^+l)!Pf(^-«)0f, 
where Ve is the Legendre polynomial, and k.h the angle between photon momentum and comoving 
wavenumber. In the limit of zero Thomson optical depth we get 

arj 

d0i k kifr 

— -Wo - — , 

d;; 3 3 

where 77 = dt/a is the conformal time. Using the fact that in the matter-dominated era the potential is 
suppresses by a f actor of 9/10 com pared to primordial value on large scales, we have 0o = -2i/'/3 
(see for example lDodelsonll2003l) and looking for a solution with time dependence of type 0f oc 



gifdrjai ^ _ -j^i Thus free streaming damps the perturbations by a factor of g^/''''*/^ - 

^-OAjdrik^ Finally we want to mention that this solution is very approximate, and in principle we 
should also take higher i modes into account. For example, on going up to ^ = 3 we get for the 
damping factor g-"-3 /t^'; ^ Also for free streaming we have super position of blackbodies, and the 



y distortion is given by dChluba & Sunvaev 



2004 Fsz = 1/2 < 0^ 1/32 < 61 >, and so the 



equivalent energy release is AE/E = 4Ysz - 1/8 < 5^ >. 
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